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Abstract. Logarithmic spirals are conjectured to be optimal escape 
paths from a half plane ocean. Assuming this, we find the rate of increase 
for both min-max and min-mean interpretations of "optimal". For the one- 
dimensional analog, which we call logarithmic coils, our min-mean solution 
differs from a widely-cited published account. 

A ship is lost in a dense fog at sea and must reach land as soon as possible. 
The captain knows that the shore is straight line, but has no information about its 
distance or direction. Equivalently, the sea is known to be a half plane, but the ship's 
location and orientation relative to the boundary is unknown. Assuming its speed is 
constant, what is the best path for the ship to follow in its search for the shore? 

The word "best" can be understood in several ways pp. We start with minimizing 
the worst-case scenario (min-max); a relevant conjecture that the family of logarith- 
mic spirals contains the minimal path remains open. Our small contribution is that of 
providing the computational details that underlie a proposition due to Baeza- Yates, 
Culberson & Rawlins [3 13 Ej • We then adopt a different sense of "best" and deter- 
mine the logarithmic spiral that minimizes the expected pathlength (min-mean), in 
which shoreline directions are assumed to be uniformly distributed. Except for the 
(admittedly large) theoretical gap regarding the optimality of logarithmic spirals, the 
calculations in this two-dimensional setting are straightforward. 

We subsequently turn to the one-dimensional analog of the search problem. The 
shore is now simply a point on a line and the spirals here are necessarily self- 
intersecting. A large computer science literature on this problem exists. The solution 
of the min-max problem was first found by Beck & Newman fS]. Their approach to 
the min-mean problem, however, suffers from the assumption of a nonuniform target 
distribution (a certain scaling property, true in the two-dimensional setting, is less 
apparent here). We give our solution, which is distinct from theirs, and hope to 
initiate discussion on this issue. 

The three-dimensional analog, for which shores are planes in space, would seem 
to be very difficult. We wonder if an appropriate extension of spiral has ever been 
examined in the past. 

^Copyright © 2005 by Steven R. Finch. All rights reserved. 
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0.1. Plancir Setting: Min-Max. Let k, > 0. Three preliminary results are: 
Lemma 1. The distance between the hne Ax + By + C — and the origin is 

\c\/VWTb^. 

Lemma 2. The equation of a hne tangent to the circle of radius R, center at the 
origin, is r — Rsec{9 — cu), where cu corresponds to the point of tangency. 

Lemma 3. The equation of a line tangent to the spiral r = e'^^ is y — e'*^sin(^) = 
m{x — e'^^ cos{9)), where 9 corresponds to the point of tangency and the slope is given 

by 

Ksin(^^) + cos(^^) 



m — 



K cos{9) — sin(6') 



Proof of Lemma 1. The unit vector (A, B) / ^/ + B"^ is normal to the line 
Ax + By + C = 0, hence the point {—CA, —CB)/{A^ + B"^) on the hne determines 
its distance from (0,0). 

Proof of Lemma 2. In rectangular coordinates, the hne is given by 

y = _Rsin(c<j) — cot(ti;) {x — i?cos(a;)) . 
In polar coordinates, therefore, wc have 

rsin(^) = Rsin{u!) — cot(a;) (rcos(^) — i?cos(a;)) 

and so 

r sin(cj) + cot(c<j) cos(ci;) , . 

R " sin(^) + cot(a;) cos(^) " ^^^^ ~ 
Proof of Lemma 3. Clearly 

dy _ dy/d9 _ (e''^ sin(^))' _ k sin(^) + cos(^) 
dx ~ dx/d9 ~ (e«^cos(6'))' ~ Kcos{e) -sm{9)' 

Theorem 4. Of all lines tangent to the spiral r = e"^, there is exactly one that is 
tangent to the circle of radius R, center at the origin. Call this line L. The tangency 
point of L with the spiral is 

^o = ^(ln(i?) + ^ln(l + «2)). 

The tangency point of L with the circle is 

uJo^9o- arccos ( , ^ ) < 6*0; 



thus L has equation r = Rsec{9 — ujq). 
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Proof of Theorem 4. Apply Lemma 1 with A = m, B = —1 and C = e^^{sm{6) — 
mcos{6)) to obtain as an expression in k, 6, m. Lemma 3 further gives m as an 
expression in k, 6. We find that + k^) = e^'^^ and hence the formula for 6q is 
true. By Lemma 2, e'^^° = i?sec(6'o — uj) and thus the formula for ojq is true. 

Think of a ship, starting at the origin and moving along the spiral. Its first contact 
point with L is at 6*0. We wish to compute its second contact point 9i. The reason 
is that, in the interval 6o < 6 < 6i, the spiral intersects all other tangent lines to 
the circle of radius R. At 9 = 9i, repetition begins so we stop there: All possible 
shorelines at distance R from the origin have at this point been found. 

A closed-form expression for 6i is not known, but it uniquely satisfies the equation 

e'^^i = i?sec(ei -cuo), < < + 27r. 

Once we have 6i for R = 1, we have it for all R via the formula 

e,{R) = -\n{R)+e,{l) 

K 

since 

uJo{R) = -ln(i?)+u;o(l), 

K 

SO 9i{R) - uJo{R) = 9i{l) - cjo(l) and thus 

_ Re'^^iW = i?sec(^i(l) - cuo(l)) = /?sec(^i(i?) - uJo{R)). 

Such scaling behavior is valuable here - we may consider R = 1 without loss of 
generality - but this property fails in Section 3. 

Lemma 5. The arclength of the spiral r = e'^^ up to O is 





J K 

— oo 

Proof of Lemma 5. From dr = Kc'^^dO, we deduce that ds'^ = r'^dO^ + dr'^ = 
e'^^^de^ + K'e^^'de^ = (1 + K')e'^'de\ 

With the assumption that i? = 1, the two-dimensional min-max problem re- 
duces to minimizing as a function of k. While an explicit for- 
mula for 9i in terms of k is unavailable, a purely numerical scheme suffices to give 
K = 0.2124695594... with arclength 13.8111351795.... The latter is consistent with the 
estimate 13.81 reported in |2]; earlier estimates 0.22325 and 13.49 from |3] arose 
when erroneously minimizing e^^^K. 

We now obtain trigonometric equations that serve to define the best spiral more 
precisely. 
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Figure 1: Two helpful pictures for the proof of Theorems 6 and 7. 
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Theorem 6. The min-max logarithmic spiral has parameter k, — tana = 0.2124695594... = 
ln(l. 2367284662...) with arclength cscasec/? = 13.8111351795..., where a, (5 satisfy 
the simultaneous equations 

1 1 271 — a — 6 cos a ^ m+o^^ 

I 1^ -,(27r— a— p) tan a 



tan a tan P cos^ a cos P 

Proof of Theorem 6. Define angles a, (3 and lengths u, v by 

^0 = « + u = = sec a, 

9i = {2% - a - p) + 00, V = e^^^ = sec p. 

Differentiating with respect to a, we obtain 

u' — secatana — utana, (1) 

v' — P' sec P tan P — P'v tan 

that is, 

v' 

P' = -cotp = v'cospcotp. (2) 



V 

From e'^'^i = g«eog«(27r-a-/3)^ follows that 

^; ^-ue^^'^-"-^)*^''" (3) 

hence 



V 



{2tt -a- P) sec^ a - (1 + tana) v + u' - (4) 



= {^—v' tan a cos /3 cot P + (27r — a — /3) sec^ a) v 
by (1) and (2). Since the objective function 

e ^ _ y (^g(^ Q, 

is minimized when 

v'csca — i;cscacota = 0, (5) 
we have the additional formula v' — v cot a. Substituting this twice into (4) yields 

cot a — —V cos P cot P + (27r — a — P) sec^ a, 

thus 

cot a + cot /3 = (27r — a — /9) sec^ a. 
Also, (3) implies immediately that 

sec/^^secae^^'^-"-^)*^'^" 

which we observe is true independent of (5). 
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0.2. Plancir Setting: Min-Mean. On the basis of Section 1, the two-dimensional 
min-mean problem clearly reduces to minimizing the average arclength 



1 ^/l 



■ Lj()+27r 



2ti k 



27rK 
2t:k 



1 - 



Ve2«« - 1, 



de 



Ve2«^ - 1, 



d9 



as a function of assuming R= 1. Define for convenience 

P) — (—2 CSC a + ln(sec a + tan a) + ln(sec (3 + tan (cot a + cot (3) , 

\E'(a, /3) = (a + /? — 27r) (sec a esc /5 + esc a sec /3) sec a, 
5(q;, /9) = sec a — cot a esc /3 + (tan a cot /3 — esc a esc /3) sec a — (cot^ a + csc^ a) sec /3. 

Theorem 7. The min-mean logarithmic spiral has parameter k, — tana = 0.3732051316.. 
ln(1.4523822387...) witi arclength 

— (ln(sec a + tan a) + ln(sec /? + tan /3) — (sec a — sec /?) cot a) esc a = 7.0321857865..., 



where a, (3 satisfy the simultaneous equations 



cos a 



cos P 



~ e 



{2-K—a—j3) tana 



Proof of Theorem 7. By the observation at the end of the proof of Theorem 6, the 
truth of the second equation is independent of the objective function. It remains to 
derive the first equation. Define 



w 



(v — u) cot q; + In + y/v"^ — 1^ +\n(u + Vu'^ — 1^ , 



then 



w' — {v' — u') cot a — {v — u) csc^ a + 



= v' cot a — u — iv — u) csc^ a + 



u 



(6) 



+ 



Titan a 



1 
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using (1). The objective function 



K 



+ ln e 



K 



+ In ( e^^'' + \/e2''^o _ i 



w CSC a 



27r 



is minimized when w' = w cot a, as with (5). Between (6) and this additional formula, 
we eliminate w' and solve for v' in terms of a, /3, m, v. Substituting the resulting 
expression for v' into (4) gives an equation involving $, \E' and S. 

0.3. Linear Setting: Min-Max. Let 7 > 1. The one- dimensional analog of the 
logarithmic spiral r = e'^^ we study here is 

x = (-7)L*J (l-(7 + l)(t- W))- 

For lack of standard phraseology, we call this a logarithmic coil. Local maximum 
points occur at (x, t) = (7^*, 2i) where i is an integer; local minimum points occur at 
(a;,t) = (-7'-\2z-l). 

Given a point X > 0, the distance 5 that the ship travels to reach X is 



2i+l 



5 = X + 2 J2 1^ = X + 



27 



2i+2 



j=-oo 



7 



where 7^* < X < 7^*"^^, that is. 



rin(X) 1 




1 


rin(X) 1 




21n(7) 




2 


ln(7) 





Given a point X < 0, the corresponding distance 6 is 



2i 

6 = -X + 2 J2 7^ 

j = -OD 

where 7^*^^ < —X < 7^*"*"^, that is, 

In(-X) _ 1' 
~ 2 



-X + 



27 



2i+l 



7 



21n(7) 



In(-X) 
ln(7) 



- 1 



Scaling as observed in Section 1 no longer works here: for the points ±1, we have 



5(1) = 1 + 



5(-l) = l + 



27 



7 — 1 ' ^ ' ' 7 — 1 

which are not easily related to S{X) and 6{—X). An analysis of S{X) and 6{—X), 
which possess sizeable jump discontinuities at 7^* and — 7^*^^, would seem to require 
different tools than before. Computer scientists traditionally normalize by |X|; see 
Figure 2 for a sample result. For more on the following theorem, see |5| l6| ITf IHj l9| ITn| 
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Figure 2: Graph of 5{X)/\X\ for 7 = 2. 



Theorem 8. The min-max logarithmic coil has parameter 7 = 2 with worst-case 
ratio 5/\X\ = 9. 

Proof of Theorem 8. For simpHcity, we examine only positive X. If X — for 
some small £ > 0, then i — k and 5/X ^ (27^ + 7 — l)/(7 — 1) as £ — > 0"^. Calculus 
gives that 7 = 2 is the critical point, which yields in turn the least maximum value 



0.4. Linear Setting: Min-Mean. The use of 5/|X| in defining the min-max coil 
in Section 3 seems fairly natural; the formulation behind a min-mean coil, however, 
requires some careful thought. Consider the integral 



whose graph appears in Figure 3. (The minimum and maximum values are 1 + 
61n(2) — 5.1588... and 1 + 12/e = 5.4145... when 7 = 2.) We wish to determine 7 for 
which /, a kind of normalized average, is minimal. Since / itself is a periodic function 
of X (although smoother than the word "minimal" can be used only loosely. 

Assuming 7^* < X < 7^*+^ and — 00 < p < i, we have 



5/\X\ = 9. 




-X 
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Figure 3: Graph of I{X) for 7 = 2. 



where 

f{i, 7, X)^{X- 7^^) + (ln(X) - 2i ln(7)) . 

Assuming — 7^^'+^ < < —7^^"^ and —00 < q < j, we have 

' /• 1 / 2729+1 \ ^ /• 1 / 272^+1 \ 



_-y2q+l 

where 



Clearly 



1 



W = ^| E /(P,7,7'^+')+/(^,7,^)- E ^(?,7,-7'^+')-^(j,7,-^) 



\p=—oo q=—oo 



and, if [ln(X)/ln(7) — 1] is even, then i = j. Upon summation, it can be proved 
that the minimum and maximum values of I{X) are, respectively, 

(7—1)^ 67—1 

The former quantity is least when 7 = 5.7041372673...; the latter quantity is least 
when 7 = 3.2232549401.... The corresponding mean ratios are 4.0089813375... and 
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4.8131558458.... These values together constitute our solution to the min-mean prob- 
lem. 

An alternative approach is due to Beck & Newman 13 IHl UHl Hlj ■ It uses a 
single random variable if, assumed to be uniformly distributed on the interval [0, 2), 
to sample different logarithmic coils with rate of increase 7. For simplicity, take 
X > 0. Then 

/ 2i+l 

E{S{X)) = X + 2 E E t'^"" I t'*^"" < ^ < 7''^'^'' 

\j=-oo 

7 - 1 V T 

= X + ^e(X7^-^) 
7 — 1 ^ / 

= X + ^/x7-^-c^/. = xfl + I±i') 
7-17 2 V M7)/ 

and this is least when 7 = 3.591121476669... = l/W{l/e), where W denotes Lam- 
bert's function. A search strategy as such is called a mixed strategy (in game theory) 
or a random strategy (in computer science). Note, however, that a uniform distribu- 
tion on H does not imply a uniform distribution on X. It is not clear to us whether 
the optimal mixed strategy (with 7 = l/W{l/e)) is necessarily preferable to our 
deterministic strategy discussed earlier. 

0.5. Acknowledgements. We grateful to Ricardo Baeza- Yates, John Shonder 
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